Experimental Lagrangian structure functions in turbulence 
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Lagrangian properties obtained from a Particle Tracking Velocimetry experiment in a turbulent 
flow at intermediate Reynolds number are presented. Accurate sampling of particle trajectories is 
essential in order to obtain the Lagrangian structure functions and to measure intermittency at 
small temporal scales. The finiteness of the measurement volume can bias the results significantly. 
We present a robust way to overcome this obstacle. Despite no fully developed inertial range we 
observe strong intermittency at the scale of dissipation. The multifractal model is only partially 
able to reproduce the results. 
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I. INTRODUCTION 

Turbulent flow still continues to puzzle. Whereas the 
study of turbulence in the Eulerian framework has gone 
through a fruitful period of discovery the main focus to- 
day is on the more subtle Lagrangian properties of fluid 
particle behavior. Global issues such as dispersion of 
pollutants, cloud dynamics, oceanic food chain dynamics 
and a variety of applications ranging from aerodynamics 
to combustion need accurate modelling and demand the 
latest knowledge on Lagrangian behavior. During the 
last 10 years we have experienced a bloom in studies of 
Lagrangian statistics of turbulence in fluid flow. Direct 
Numerical Simulation (DNS) studies [H, i, i, i, i, [| 
and laboratory experiments (inainl y P article Tracking 
Velocimetry (PTV)) 0, i, & IIl|,M El, Q P have 
played an important role in revealing the physics gov- 
erning the behavior of particles at the smallest scales in 
a turbulent flow. Along with this, the theoretical un- 
derstanding has improved: the multifractal model p^ . 
originally introduced in the Eulerian framework, has now 
turned into a promising phenomenological model for La- 
grangian observations |17| . 

In this contribution we present a Lagrangian analy- 
sis of small scale statistical behavior through higher or- 
der velocity structure functions. From a PTV experi- 
ment we obtain particle trajectories and from these we 
construct structure functions of velocity along the tra- 
jectories. This exercise is common in the field of turbu- 
lence and the above mentioned references all have the La- 
grangian structure functions as the starting point. Stud- 
ies of the smallest time scales of the flow has revealed 
intermittent behavior. The first signs were observed 
with DNS. Only within the last couple of years has it 
been possible to measure Lagrangian intermittency in 
a physical flow with PTV and hence quantitatively de- 
scribe the extreme statistics present in the Lagrangian 



data |,i[Tl,[Tl,[Tl. 

The joint work by the International Collaboration of 
Turbulence presented in Arneodo et al. [l3| showed that 
the big picture is the same, whether you use DNS or 
PTV. DNS and physical flows does, however, have both 
quantitative and qualitative differences. DNS has the 
disadvantage that at present, due to limited computer 
power, the largest structures in the flow can only be sim- 
ulated a few times, and hence the statistics becomes very 
poor on large scales. As we will argue in this paper, this 
could have an influence on even the smallest structures 
in the flow. In contrast we can do many independent 
realizations of the flow with PTV and make long runs 
so that all scales are well resolved and statistically well 
represented in the many ensembles. Unfortunately, mea- 
suring in a finite volume of the flow can bias the results 
significantly: fast particles will leave the volume early 
and hence statistics for long times are primarily based on 
slow particles. These differences along with a few more 
should not be neglected when analyzing data since wrong 
conclusions could then easily be made. These issues are 
the major motivation behind this paper. We will present 
a thoroughly way through the jungle of random errors 
and bias and try to quantify the importance of each. 

The paper is structured as follows: the technique of 
PTV is explained and the flow is characterized in Sec- 
tion |TT1 In Section IlIII we present a robust way to quan- 
tify at which scale the finite volume bias sets in. We 
calculate the final structure functions in Section IIVI and 
finally make a comparison with the multifractal model. 

Central to the paper are the Lagrangian structure func- 
tions S'p(r) of order p. With the velocity increment along 
a trajectory 



(5,v(t) = v(i + T) - v(t), 
we define a structure function of order p as 
5p(r) = (|e.<5.v(T)n, 



(1) 



(2) 
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with (•) denoting ensemble averaging as well as averaging 
over all possible directions of the axes, i.e. the random 
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unit vector e: 

Sp{r) = (|,5.v(T)ncos(5.v(r),e)r) (3) 
= (|<5.v(T)n(|cos(<5.v(r),e)r) 



where we have used 

/ d6lsin6'|cosP6l| = — (4) 
Jo ' 1+p 

We thus present >S'p(r) in a way independent of any 
particular choice of coordinate system since we have in- 
cluded all possible rotations around the center of some 
spherical volume. With this definition the ensemble is 
close to isotropic. For the flow studied in this paper 
isotropy is actually only strictly present in the center of 
the tank. The use of this ensemble does have several ad- 
vantages. We can compare with other experiments and 
DNS regardless of any degree of anisotropy and hence 
move closer towards a general understanding of any uni- 
versal behavior regardless of anisotropy; deviations from 
theories such as the multifractal model will not be ex- 
plained by the presence of anisotropy. Since isotropy is a 
key element in the foundation of the multifractal model 
any such analysis must therefore refer to isotropic be- 
havior and hence isotropic ensembles of data. Attempts 
to develop a framework for anisotropy in turbulence has, 
however, been done [13, 

The time scales relevant for particle motion range from 
the viscous scale (the Kolmogorov time scale) to the 
integral time scale T^. The Reynolds number, measuring 
the strength of the turbulence, scales like Re\ ~ Tl/t,,. 
Two-time particle statistics, where the time lag t — ti — 
t2 is less than the integral time scale, but larger than the 
Kolmogorov time scale r^, are said to be in the inertial 
subrange. From an experimentalist's point of view this 
means that Lagrangian inertial subrange scaling is very 
difficult to obtain compared to Eulerian statistics where 
the size of the inertial subrange grows as Re^^^ . 

K41 similarity scaling predicts that for time lags in 
the inertial range S'p(t) ~ e^'/^r^/^, where e is the mean 
kinetic energy dissipation. The lack of similarity intro- 
duces corrections to K41 similarity scaling. These are 
due to intermittent events. The multifractal model [iGj 
is today the most used model of intermittency. The lack 
of an inertial subrange of Splr) observed in experiments 
and DNS indicates that the Reynolds number is not the 
crucial factor for observing intermittency and data from 
a range of Reynolds numbers have also been observed to 
follow each other closely [l7| . 

II. PARTICLE TRACKING VELOCIMETRY 
A. Experimental Setup 

We have performed a Particle Tracking Velocimetry 
(PTV) experiment in an intermediate Reynolds number 



turbulent flow. PTV is an experimental method suitable 
for obtaining Lagrangian statistics in turbulent flows. 
Lagrangian trajectories of fluid particles in water are ob- 
tained by tracking neutrally buoyant particles in space 
and time. The flow is generated by eight rotating pro- 
pellers, which change their rotational direction in fixed 
intervals in order to suppress a mean flow, placed in the 
corners of a tank with dimensions 32 x 32 x SOcm'^ (see 
FigllJ. The data acquisition system consists of four coni- 




FIG. 1: Experimental setup 



mercial CCD cameras with a maximum frame rate of 
50Hz at 1000 x 1000 pixels. The measurement volume 
covers roughly lOOOcm'^. We use polystyrene particles 
with size ~ 400/im and density very close to Ig/cm"^. We 
record 0(800) particles at each time step with an accu- 
racy in the estimation of the particle position of 0.05 pix- 
els corresponding to a standard deviation ajuter ~ 10/im. 
The particles are illuminated by a 250W flash lamp. 

The Stokes number St = Tn/rri measures the ratio be- 
tween the relaxation time, r^, of particle motion relative 
to the fluid and the Kolmogorov time scale, r^. Here 
tr — {^/^^){Pp/ Pf){dp/j^) where p/ is the density of the 
fluid and pp and dp are the density and the size of the 
particles respectively. We get St = 0.01 and thus much 
less than one. The particles can therefore be treated as 
passive tracers in the flow. 

The mathematical algorithms for translating two di- 
mensional image coordinates from the four camera chips 
into a full set of three dimensional trajectories in time 
involve several crucial steps: fitting 2d gaussian profiles 
to the 2d images, stereo matching (line of sight crossings) 
with a two media (water-air) optical model and construc- 
tion of 3d trajectories in time by using a kinematic prin- 
ciple of minimum change in acceleration [22l . |23| | . 

If a particle can not be observed from at least three 
cameras the linking ends. Most of the the time this hap- 
pens because the particles shade for each other. The 
higher the seeding density of particles the higher the risk 
of shadowing. Since the seeding is relatively high we ob- 
tain shorter tracks than we would have in the case of 
only a few particles in the tank. Since the particles most 
often only disappear from one or two cameras for a few 
time steps, a new track starts when the particle is again 
in view from at least three cameras. The track is hence 
broken into smaller segments. Through kinematic pre- 
diction we are able to connect the broken track segments 
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into longer tracks. We use the method of Xu (2j|. The 
result is satisfactory with a substantial increase in the 
mean length of the tracks. 

The flow characteristics are presented in TablelH With 
= 0.08s and a recording frequency at 50Hz the tem- 
poral resolution is ^ 4frames/r^. The mean flow is ax- 
isymmetric with a significant vertical straining on the 
largest scales and no significant differences from the flow 
reported in [isl . where properties of the mean flow 
can be found. 

We choose a coordinate system centered approximately 
in the center of the tank where the velocity standard de- 
viation (T„ has a global minimum. The radial distribution 
of particles in the measuring volume is presented by the 
stars in Figure [21 We see that for distances from the cen- 
ter less than 50mm, represented by the vertical dashed 
line, the particles are uniformly distributed. We therefore 
choose a ball, B, with same the center and a radius of 
50mm, as the volume for all future studies in this paper. 
In addition wc see that both the start and end position 
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FIG. 2: Pdf p(r) of the radial position of particles in the measuring 
volume (* symbol). The symbols □ and o represent respectively 
particles starting and ending a trajectory. The black solid line has a 
slope equal to two and therefore represent the uniform distribution. 
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FIG. 3: Top: PDF histogram of number of particles within ran- 
domly positioned balls of radius 10 mm (left) and 30 mm {right). 
The dots represent a Poisson distribution. Bottom: Variation of the 
ratio ({N'^) — (N)^)/{N) as function of ball radius. The horizontal 
line represent the Poisson value of unity. 
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TABLE I: Turbulence characteristics: e is the mean kinetic energy 
dissipation, rj = (y^ je)^!'^ is the Kolmogorov length scale with the 
kinematic viscosity v = 0.89mm^/s of water, t,, = (u/e)^^^ is 
the Kolmogorov time scale and cr^ = ^ (fu^ + o"^^ + ct^^ ) is the 
standard deviation of velocity. The integral length scale is defined 
as L = a^/e while is the eddy turnover time = L/au- The 
Reynolds number is defined as Rex = with the Taylor micro 



scale A : 



of trajectories are also uniformly distributed within the 
ball, B. This means that tracking failure is independent 
of position. Trajectories may move in and out of B but 
only positions inside the sub-volume are considered in 
the data analysis. 

To see whether the particles are statistically indepen- 
dent we check the data against a Poisson distribution. 
For every tenth frame we place 100 balls of varying ra- 
dius randomly within the flow and count the number of 
particles inside. The result is presented in Figure [3l The 
top figures show two examples; one for particles inside 
balls of radius 10 mm and one where the ball radius is 
30 mm. From the bottom figure we conclude that the 
particles obey a Poisson statistic within 8% which is a 
bit better than the PTV experiment presented in [25j . 
This means that the particles either cluster or repel each 
other. 



The database of trajectories is compiled from 73 runs 
performed under identical conditions. Each run consists 
of 10000 consecutive frames. After the recording of a run 
the system was paused for three minutes before a new 
run was recorded. We therefore consider the 73 individ- 
ual runs to be statistically independent. Throughout the 
paper error bars will therefore be calculated as statistical 
standard errors of the mean. 



B. Binomial filtering 

Even though we know the position error epos to be 
very small, we choose to filter the data. We choose a bi- 
nomial filter which has the convenient property of com- 
pact support. Using a binomial filter instead of a conven- 
tional Gaussian, does not seem to have a large effect (not 
shown). The weights Wk in a binomial filter of length N 
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is given by 

Wfe = 2 



1-^'^^ k = 0,l,...,N-l 



(5) 



The width of the filter is af^uer = At^/N ~ 1/2. We 
apply the binomial filter on the position measurements 
treating each dimension separately. The velocity and ac- 
celeration are then simply given by finite differences 



v^{t + At/2) 



Xi{t + At) - Xi{t) 
Xt 

x,{t + At) ~ 2xi{t) +Xi{t- At) 
Ai2 ' 



(6) 



where ~ denotes filtered quantities. 

We inspect filters with a length from iV = 1 (unfil- 
tered) to TV = 20 and look at the standard deviation and 
flatness of velocity and acceleration as functions of the 
filter standard deviation relative to the Kolmogorov time 
scale Tr^. These are displayed in figured] While velocity 
seems to be unaltered by the filtering the acceleration is 
very dependent on filter length. This is a general problem 
with measurements of accelerations. Both the standard 
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FIG. 4: left: ctu {*) and aa (•) as functions of filter width a-futer 
divided by Kolmogorov time scale r^. right: Flatness Fu (*) and 
Fa {•). The dotted vertical lines denote a filter of length N = A 
corresponding to 0.23t,j. 



deviation and flatness are functions of the filter width. 
The amount of filtering is a trade off between eliminat- 
ing noise and eliminating the real signal. The figure does 
not give a clear indication of which filter to use. We 
therefore also look at the acceleration pdf. Since this is 
the limiting PDF for velocity increments i5v(r) it might 
seem important for the outcome of results on Lagrangian 
structure functions. In figure [5] we show the pdf of accel- 
eration as a function of the filter length. From = 4 the 
shape of the pdf becomes more or less constant with N. 
The fat tails of the pdf based on unfiltered data {N = 1) 
are likely to be due to noise and perhaps bad connection 
of tracks. Filtering removes the noisy tails and makes 
it possible to estimate moments up to (a^). The critical 
filter length for which convergence is achieved is iV = 4. 

= 4 corresponds to a filter width of 0.23r^, i.e. well be- 
low the Kolmogorov scale. It should be emphasized that 
choosing = 4 is still somewhat arbitrary since there is 
no rigorous way to determine the optimal filter. Unless 
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FIG. 5: Acceleration pdf p(ay) for the component along the axis 
of forcing in the tank. The pdf is normalized with the standard 
deviation of acceleration for the same component. The differ- 
ent colors represent the different filter lengths ranging from 1 to 
20 following the direction of the arrow. The two black curves are 
the gaussian distribution and a stretched exponential fitted to the 
curve for N = A (only partly visible as the thin black curve below 
the fitted stretched exponential). 



otherwise stated N = A has been used in the remainder 
of the paper. 

We have fitted a stretched exponential of the form 
p{a) = 7V/exp(-aV((l -I- \a(3f/af\)'>f)aj) to the pdf of 
scaled acceleration for N — A. With f3f — 0.65, cr/ — 0.55 
and 7/ — 1.40 the fit is excellent. The same functional 
fit for p{a) was used by Porta et al. Q with slightly dif- 
ferent fit parameters. We can also give an estimate for 
the dimensionless constant qq in the Heisenberg-Yaglom 
relation (aiaj) = aQe^^'^v~^l'^8ij. In principle ao is not 
a constant but a function of Reynolds number. It has 
been well studied in the literature and it is closely con- 
nected to Lagrangian stochastic models. Sawford et al. 
[1^ fits functional forms of ao as function of Reynolds 
number to data obtained from DNS and high Reynolds 
number PTV. The values obtained from the PTV data 
are approximately 30% larger than for the corresponding 
DNS data. In addition strong anisotropy is observed in 
the PTV data. We get ao — 4.5 ± 0.3 for the component 
averaged acceleration variance while only oq = 3.8 ± 0.2 
for the component along the axis of forcing. For the axial 
component f2(| suggest oq = 6.5(1 + 134i?e^^)^^ giving 
ao = 3.3 for the present Reynolds number. The con- 
struction of our apparatus does, however, not give us the 
opportunity to investigate a large span in Reynolds num- 
bers which would be necessary to verify functional forms 
of ao. 

As already mentioned the filtering has a large impact 
on acceleration. Besides that, there is a slight chance that 
we might underestimate accelerations since variations at 
the smallest temporal scales sirnply cannot be resolved. 
In the PTV data presented in [2y| the ratio between 

and sampling frequency is ~ 23 frames/r^ (for Re\ ~ 
870) compared to only 4 frames/rr^ in our experiment. 
Their data is, however, much noisier, so that their chosen 
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filter width of u f uteri '<'r\ = 0.15 is close to our of 0.23 for 
iV = 4. 

The motivation behind filtering the trajectories was to 
eliminate noise, i.e. the error, epos, associated with deter- 
mining the position of a particle. Even though epos might 
be both random, unbiased and uncorrelated it still con- 
tributes to, for example, the Lagrangian structure func- 
tions. To see this we assume that Xi(t) = Xi{t) + epos{t). 
Xi(t) is a measured component of position on a trajectory 
while Xi (t) is the true position. Furthermore we assume 



posit)) = and {epositi)eposit2)) 



a 



jitter 



if h - t2 



and zero otherwise. 

The correction to the structure functions is now 
straight forward to calculate. Since a different number 
of particle positions are involved in the calculation as a 
function of time lag, r, the correction becomes a function 
of the time lag, r, itself. With filter = 4 we get for the 
second order measured structure function 5™(t) 



where cr?„ 



STir)^S2{r)+a. 
is given by 



2 

error 



(r) 



(7) 
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(8) 



We find that the correction to the measured higher order 
structure functions S'™(r) is given by 

Sl^ir) = 5p(T) + Apal,,,{T)S:;_^{T), (9) 

with 5*0 = 1 for all time lags r. We have omitted terms 
of higher order in a err or- -^p is given by 



(10) 



With a jitter ^ 10/im the corrections cr^^-rori''') ^^'^ of or- 
der 10~^mm^/s^. 



III. BIAS OF LAGRANGIAN STATISTICS 



is dissipated and hence we have ^^(v^) ~ — e [271. The 
equation can be derived directly from the Navier-Stokes 
equation by assuming global homogeneity [20| which, as 
already mentioned, is only approximately true for this 
experiment. 

In DNS the random forcing occurs in fourier space 
and is hence globally homogeneous. We therefore have 
d{v'^)/dt = and consequently Lagrangian stationarity. 
However, most physical flows encountered in nature, as 
for example the atmospheric boundary layer, will seen in 
a finite volume be Lagrangian non-stationary. 

Particles with fast velocity tend to leave the measure- 
ment volume after only a short amount of time. This 
means that those particles that stays in the volume for 
long times often are those with the smallest velocity. 
The effect is a bias for long times towards slow parti- 
cles. It should be emphasized that this is a systematic 
error whereas the Lagrangian non-stationarity is a gen- 
uine property of the flow. 

Exactly how one should compensate for the systematic 
error is an open question. In this paper we will build on 
the ideas first presented by Ott and Mann [131 . Here 
the Lagrangian structure functions S'p(t) are expressed 
through the mean Greens function G(r, r) as 



Spir) 



Sp{T\v)G(r,T)d' 



(11) 



where Sp{t\v) is the conditional Lagrangian structure 
function of order p defined as the mean of \5v\p condi- 
tioned on the distance travelled r after a time lag r. The 
mean Greens function for one-particle diffusion G(r, r) 
is the probability density of getting r after a time lag r. 
The relation is evidently exact. However, when the field 
of view is limited the integration on the right hand side 
is truncated which leads to systematic errors. It is trun- 
cated by a filter W{r) that expresses the probability that 
a point X2 separated a distance r from another point xi 
lies inside a ball with radius R. Here Xi is chosen ran- 
domly inside the ball. Assuming homogeneity W(r) is a 
purely geometric factor of the distance |xi — X2I = r. It 
is given by 



W{r) = 



2R' 



AR) 



for r < 2R, 
for r > 2R. 



(12) 



A practical property of the present experiment is the 
stationarity of Eulerian velocity statistics. The La- 
grangian statistics are, on the other hand, not stationary 
in the measurement volume, B. 

A particle which enters B will loose kinetic energy dur- 
ing its travel inside B. This reflects the non-uniform forc- 
ing in space in our experiment. On average the particles 
gain kinetic energy close to the propellers located outside 
B. During their subsequent motion the particles lose ki- 
netic energy until they again come close to the propellers 
which are constantly spinning. Thus there is a flux of ki- 
netic energy into B. Inside the volume the kinetic energy 



If we neglect the finite measurement volume and just 
average the velocity differences that we have actually 
measured, the experimental structure function becomes 



5. 



p^meas 



(r) 



J W{r)Sp{T\r)Anr^G{r, T)dr 
/iy(r)47rr2G(r, r)dr ' 



(13) 



lOtt and MannI proposed an improved method where 
W{r) is removed from eqn. 1131 each pair is binned with 
weight 1/W{r). Since W{2R) — 0, pairs with separa- 
tions very close to 2R should be disregarded since they 
would otherwise make the compensation explode. We 
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therefore limit separations to 2R — S. In the present case 
2R — 100mm and d — 5mm. 

Including compensation we can now write eqn. [13] as 




20 40 60 80 100 

r [mm] 



FIG. 6: Mean Greens function G(r, t) as a function of r at three 
different time lags r. Increasing towards the right, time lags are 
T = {5, 13, 35}t,, with (solid lines) and without (dashed lines) com- 
pensation through W{r). 



We will use this framework to estimate an upper time 
lag below which Lagrangian structure functions are not 
biased by the finite measurement volume. To put it sim- 
ple we need to figure out whether or not 2R — d is enough 
to cover the support of the integrands in eqn. (M) 

In Fig. [5] we show the mean Greens function G{r, r) as 
a function of r. For the small time lag r = 5r^ the curves 
almost collapse indicating that finite volume effects are 
almost vanishing. Increasing r to 13t,, we see that finite 
volume effects are significant as the curves no longer col- 
lapse. Both curves, however, go to zero for large values of 
r, which means that the integrand in eqn. 1141 converges. 
This means that for this specific time lag we can calculate 
finite volume unbiased structure functions if we include 
compensation. For the largest time lag, r = 35r^ we see 
that the compensated G(r, r) no longer converges. This 
means that even the compensation fails. 

We now assume G(r, r) to be Gaussian and self simi- 
lar, that is G(r,r) = exp(-(r/a,{r))V2)/(27r3/2a,(r)3). 
Furthermore we assume that Sp{T\r) ~ for large r 
and large r. Figure [7] shows log-log plots of Sp{T\r) 
for p = 2,4,6,8. We observe an approximately agree- 
ment with Sp{T\r) ^ r^. Ott and Mann [27| calculated 
S2 for a gaussian displacement process and found that 
S2(T\r) ^ for large r. The implications for Sp is ev- 
ident: only the slightest deviation from zero in G(r,T) 
for large r will cause the integrand in the numerator of 
eqn. [14] not to converge. With the made assumptions we 
can calculate the relative error on Sp from the finite vol- 
ume. For the compensated structure function the relative 
error given by 1 — Sp^meas/ Sp. Sp,meas can be calculated 




- / 1 -z r / 1 
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FIG. 7; log Sp{T\r) / Sp{t) as a function of logr/axitau) for 
r e (15; 25)r^. The four panels represents p = 2, 4, 6, 8 with 
fits r*" for large r / Gxitau). 




0-,/R 

FIG. 8: Relative error estimates for Sp for p = 2,4,6 and 8. 
Dashed lines are without compensation while solid lines are 
with compensation 

from eqn. 1141 while Sp is given by eqn. II II integrating all 
the way to infinity. Likewise we can calculate the error 
without compensation. In this case Spiraeas can be calcu- 
lated from eqn. [T3] The results are presented in Figure [5] 
for p=2,4,6 and 8. The solid lines represent the compen- 
sated structure functions while the dashed lines represent 
the uncompensated. The error is a function of the non 
dimensional variable Ux/ R- 

It is evident that without compensation we get large, 
systematic errors even for quite small values of Ux/R 
while the compensation works up to a point where the 
upper limit of integration (2R — 5) is felt, and the com- 
pensation rapidly deteriorates as we move beyond this 
limit. The upper limit on ax{T)/R defines a critical time 
lag Tcrit where measurements for r > Tcrit exhibit a 
serious, systematic error. For p = 2 the critical limit 
is (Tx{Tcrit)/R ~ 0.5 while it is ax{Tcrit)/R ~ 0.4 for 
p — S. Actually these estimates are optimistic because 
G in practice tends to have fatter tails than a Gaussian. 
This is shown in Figure [5] 
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FIG. 9: G{r,T) for r = lOr,,. The full line is a Gaussian. 



In order to get an estimate of Tcrit we set axir) auT 
leading to 



^ ^ = 15-'/'C,j R, (15) 

with C2 --0.6, C4 -0.53, Ce -0.49 and Cs -0.45 and 
the integral scale L = a^/e. L was chosen because it is 
defined by the geometry of the apparatus used indepen- 
dent of R\. From (jlSp we see that Tcrit /Tr/ increases with 
R\, thus making it easier to measure at high Reynolds 
numbers. In other respects, such as the demand on frame 
speed, it of course gets harder. If we wish to study large 
time scales, the ratio Tcrit,p/TL, where Tl — cr^/e is the 
Lagrangian integral time scale, could be more relevant 
and we can note that 



IV. INERTIAL RANGE SCALING 
A. Presentation of data 

The linear dependence of Re\ on Tl / Trj implies that a 
very high Reynolds number is needed in order to obtain 
a clear Lagrangian inertial range. Yeung \t\ concluded, 
based on extrapolations from Eulerian fields in DNS, that 
at least Rex — 600 — 700 is needed. Experimental flows 
at Rex = 1000 [18^ and Rex 815 [2I do, however, 
not show a very pronounced range with a linear regime 
in the second-order structure function, 5'2(t) = CqSt, 
and one could speculate if such a range exists at all. Co 
plays a crucial role in stochastic models [30] and has been 
shown to reflect anisotropy in the large-scale forcing [2^ . 
In Figure [TO] we present results of Cq for the isotropic 
ensemble as well as for the three directions. According 
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FIG. 10: Co : The circles represent the two horizontal directions 
while the squares represents the axial direction. The stars are the 
mean (eqn. |4)l. The arrow indicates the maximum Co = 4.46 at 
T = 2.96r„. 



In other words, finite size effects are not affected by the 
Reynolds number at large time scales. 

Using eqn. llSl with the present data we find Tcrit.2/Tri = 

18, TcritA/T,-, — 16, TcritfilTr, — 14 and Tcrit,S,lTr^ = 13. 

This is when a 5% error is accepted. Inspection of the 
data shows that at these values the (compensated) inte- 
grands in eqn. [14] are indeed just covered within 2R — 6. 
Without compensation we find very small critical limits: 
Tcrita/Tv = 3.3, Tcrit,i/T^ = 2.0, Tcritfi/r,, = 1.6 and 

Tcrit.s/Tri = 1.2. 

Other studies have also looked at the bias effect of La- 
grangian statistics |18|, i25„ The criterion in eqn. [15] 
and eqn. [TO] are the strictest yet presented in the liter- 
ature. The most important lesson is, however, not the 
limits suggested by the equations themselves, but the 
compensation: without this, even small time lag statis- 
tics are heavily biased as seen in Figure [5] 



to ^2 (t) , Co should be determined from a plateau in the 
inertial range. The inertial range is almost vanishing 
in our experiment. The ★ symbols are calculated from 
eqn. [4] The maximum is Co = 4.46 at a time lag r/r^ — 
2.96 and therefore mainly associated with small scales. 
The o symbols represent the horizontal directions (Cq = 
4.88 and Cq — 4.69) while the □ symbols represent the 
axial direction (Cq — 3.83). A rough estimate of the 
error on Cq is 0.3, originating from a 10% error in the 
determination of the kinetic energy dissipation e. The 
statistical error is essentially zero. 

It is interesting to see that the slight anisotropy in 
the forcing is manifested all the way down to r — r^. 
The propellers forcing the flow are attached to four rods 
placed in the corners of the tank. The reason for the 
horizontal components being different is probably small 
differences in the vertical placement of the propellers on 
the rods. The lack of small-scale isotropy in the current 
experiment should not necessarily be taken as a failure 
of Kolmogorov's hypothesis of local isotropy. For such a 
statement the Reynolds number in our experiment is sim- 
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ply not high enough. Other experiments at much higher 
Reynolds number do, however, all indicate that the large 
scale inhomogeneities are also present at smaller scales 
although with smaller amplitude [1^, [3l|, Is^] ■ These find- 
ings are striking and may suggest that the hypothesis of 
local isotropy and the concept of locality should be re- 
viewed [sst - 

A theory that demands isotropy, as is the case of most 
K41-like predictions can literarily not be falsified, since 
the perfect experiment with isotropic forcing and hence 
isotropy on the smallest scales can not be constructed. 
This was one of the motivations behind the construction 
of the isotropic ensemble given in eqn. [4] 

Alternatively one can calculate Co from the velocity 
spectrum. Arguments put forward by Lien and D'Asaro 
[34| state the inertial range scaling is easier to obtain in 
Fourier space through the velocity spectrum. However, 
no difference was observed in the present data set through 
such an analysis (not shown). 

We now look at the higher order structure functions. 
We want to quantify the degree of intermittency through 
anomalous scaling exponents. The structure functions 
S'p(r) for p = 2,4,6,8 are displayed in Figure [TT] It 
should be remembered that the data are heavily influ- 
enced by finite volume bias for time lags r > Tcrit- The 
most important conclusion to state from the plot is the 
evident lack of power law behavior and hence a KAl scal- 
ing regime. 

lO'f 1 
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FIG. 11: Lagrangian structure functions Spi^r) as a function of 
time lag t. p is increasing upwards with p = 2, 4, 6, 8. The different 
structure functions have been shifted vertically for clarity. 

Motivated by the lack of a clear inertial range in ac- 
cessible turbulent data, Benzi et al. 35] introduced Ex- 
tended Self Similarity (ESS). Even though it was origi- 
nally applied to Eulerian data it in can easily be adapted 
to Lagrangian. Instead of plotting S'p(r) against time 
lag T, Spir) is plotted against the structure function not 
affected by intermittency for the corresponding time lag. 

In more general terms we define in the ESS context the 
anomalous scaling exponents Cp(''") through 

^^(^)=dlog[^.(.)]- 



It has obvious advantages compared to conventional ad 
hoc power law fitting procedures |[i3]. It can, however, 
be difficult to quantify from experimental data due to the 
derivative which is very sensitive to noise. 




0.1 0.2 0.5 1.0 2.0 5.0 10.0 20.0 
t/t. 



FIG. 12: Cpi''')- Top: p = A, middle: p = 6 and bottom: p = 8 
(solid curves). The shaded areas represent the the multifractal 
predictions: the upper boundary is based on longitudinal Eulerian 
structure functions while the lower boundary is based on transverse 
Eulerian structure functions. The thick black vertical lines denote 
the critical time lags according to Section IlIII We find the large 
time lag saturation levels to be 1.53, 1.76, 1.93 for p = 4, 6, 8. The 
error bars refer to statistical errors. 



Cp(t) for p — 4,6,8 are plotted in Figure [T2l Two 
general trends are observed in all three figures. 

A dip around 2t^ < t < 5r^ is observed. Hereafter 
a plateau is reached. For increasing order p the dip be- 
comes larger and Cp saturates at a higher level in agree- 
ment with the findings and speculations by Arneodo et al. 
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[17|. The saturation levels are displayed with horizontal 
lines at 1.53, 1.76 and 1.93 respectively. At time lags 
around Tcrit we see that the error bars grow significantly 
at least for p = 8. With increasing p it even happens 
before Tcrit as also suggested in Section Hill and makes a 
plateau hard to observe. 

The dip has been associated with vortex trapping [sH] 
where particles are trapped in strong vortices with time 
scale close to r^. An example of such a particle is dis- 
played in Fig. [ini At around t ~ 50t^ the particle expe- 
riences extreme accelerations (around 10 times the rms 
value) and seems to be caught in a vortex-like structure. 

We thus support the findings already reported in Porta 
et al. Q, Biferale et al. |3Q] in studies of much more vivid 
flow. However, we demonstrate that the value of the 
Reynolds number is not necessarily crucial in order to ob- 
serve characteristic turbulence features in the Lagrangian 
frame as is also known from low Reynolds number DNS. 
Differences in DNS and physical flows are, however, no- 
table and careful analysis of both is therefore important 
in order to obtain a complete picture. 

It is also striking that the maximum of S2 (t) /t is at 
T = 2.96r^ which is very close to the position of the 
viscous dip of C4 (''')■ Whether this has any significance 
is not at all clear. If so, it is a challenge for stochastic 
models, which do not take the viscous dip into account. 



z [mm 
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FIG. 13: Top Sample particle. The color denotes the magnitude 
of the acceleration. Bottom Particle acceleration in units of the 
standard deviation aa for all three components. 



B. Multifractal model 

In Fig. [12] we have also plotted the multifractal predic- 
tion. The multifractal model is developed in the Eu- 
lerian frame by Parisi and Frisch [s^ to characterize 
the spatial structure of dissipation in turbulence it was 
later adapted in the Lagrangian frame by Borgas [ssj . 
It gives a translation between the two frames and as 
such works as a brid ge b etween them. Work presented 
in d, [1^ m, m, [39I.I40I has shed light on the issue of 
multifractals in the Lagrangian frame through a number 
of high Reynolds number experiments and DNS which 
were well captured by the theory through ESS scaling 
relationships for time lags r > . With the extension by 
Chevillard et al. [s^ the multifractal framework is now 
capable of taking into account also the vortex trapping 
behavior taking place at time lags r ~ r^. 

In [l3| it was shown for C4(''') how the multifractal 
model matches results from experimental and numerical 
data on all time scales independently of Reynolds num- 
ber. The data set used in this paper is the one denoted 
EXP 1 in that paper. Here we show that for p = 6 and 
p = 8 the multifractal model seems somewhat less per- 
fect. 

In the multifractal model the flow is assumed to possess 



a range of scaling exponents h„ 



with a certain 



probability so that the velocity difference by separation 



is Sr 



For each scaling exponent h there is a 



fractal set with a /i-dependent dimension D{h). The em- 
bedding dimension is three (r S TZ^) and hence D{h) < 3 
for all h. The probability Phir) of having an exponent 
h at separation r is therefore proportional to 3 — D{h). 
From dimensional arguments the Eulerian velocity fluc- 
tuation SrU is related to a Lagrangian time lag t. That is 
T ~ r/SrU. Recently focus on this relation has cast seri- 
ous doubt on its usage [U, [H, El] . In [43j] it is shown how 
the relation r ^ r/drU is the limiting case of something 
more general. In three dimensional turbulence they show 
that Lagrangian statistics is much more influenced by 
Eulerian integral and dissipation scales than the simple 
picture suggest, where the time lag r is only associated 
with eddies of size r. This is also what we already know 
from Figure [S] and from eqn. II II where contributions from 
all scales is included in the integral. The conclusion must 
be that good statistics of all Eulerian scales are necessary 
in order to calculate Lagrangian structure function - even 
at small time lags. At present DNS does not resolve the 
statistics of the largest scales sufficiently. PTV experi- 
ments do not have this problem. 

Following [l3| closely we can calculate the multifractal 
prediction for p — A,6, 8. We choose the same model 
constants and functional form of D{h) as in [l3| since 
these were shown to fit a large number of experiments 
and DNS simulations for p = A. 

We again take a look at Figure [121 The result in [l3| 
(p = 4) is reproduced: Following the multifractal longi- 
tudinal curve (upper) in the dip quite close, C4(''') satu- 
rates at a value close to the multifractal transverse curve 
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(lower). For p = 6 and p = 8 the multifractal predicted 
curves does not fit the data. First, the plateau in the 
data does not reach the level of the multifractal predic- 
tion before finite volume bias effects become important. 
This bias has the effect of lowering the value of Cp (''')• 
Second, for both moments the dip is shifted towards the 
left. The dip minimum does, however, seem to match the 
data. In the multifractal model there is a free parameter 
included in the model definition of the dissipative time 
scale. We have set this constant, t scale equal to 7. Its 
only job is to scale the t/t^ axis. The same constant 
value was used in [l3] in order to fit the model predic- 
tion to data for p = 4. Since tscaie must be independent 
of p, we could equally well have fitted tscaie to the data 
for p — 6 or p — 8. If we had done so we would encounter 
a bad match for p = A and p = 8 or p = 4 and p = 6 
respectively. 

Why the multifractal model fail to predict the differ- 
ent moments of Cp(''") is a relevant question. We could 
speculate that the temporal resolution in our experiment 
is not high enough to resolve the smallest scale. In order 
to investigate this we look at the importance of filtering. 
We picture the increase of filter length as a way of de- 
creasing the temporal resolution and calculate Cp{t)) as 
a function of filter length, N. For < 10 the resuls are 
presented in Figure [HI With increasing filter length we 
see that the dip region is depleted and the position of 
the minimum on the r axis is shifted to the right (the 
arrows indicate the position of the local minimum for 
N = 1 and N = 10). The shift in position is decreasing 
with increasing p and hardly observable for p = 8. The 
saturation level is more or less constant with N. Both 
the quantitative and qualitative shape of the dip region 
of Cp(''') are therefore quite sensitive to the filter length 
TV. How does this relate to the multifractal model? In 
Figure [T^ we see that the dip position of the multifractal 
model is almost constant with p, while the data show a 
shift towards smaller time lags when p is increased. On 
the other hand, we observe in Figure[T4]that for p = 4 the 
shift is towards larger time lags whereas it is almost non- 
existing for p = 8. That is, for a higher resolution, here 
pictured as a lower filter, the multifractal model would 
work even worse. 
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FIG. 14: Cpi''')- Top: p = 4, middle: p = 6 and bottom: p = 
8 (solid curves). The colors denote different filter A'^ (increasing 
upwards from A'' = 1 (red)). Arrows indicate minimum position 
for A' = 1 (left arrow) and A^ = 10 (right arrow). 



V. DISCUSSION AND CONCLUSIONS 

Analyzing Lagrangian data uncritically can lead to 
substantial biases. Building on the ideas first presented 
by Ott and Mann [23] we have found a robust way to 
compensate for the finite volume bias. We show that the 
time lags for which finite volume bias can be neglected 
are limited. We also saw that the finite volume effects 
increase for increasing order p of the structure functions. 
Observing extreme statistics might therefore be very dif- 
ficult. We are convinced that the present study and its 
consequences should be kept in mind when designing fu- 
ture experiments for measuring Lagrangian statistics: in 



a high Reynolds number flow the separation between the 
integral scale and the Kolmogorov length is very large. 
In order to follow particles and measure the Lagrangian 
structure functions for large time lags without finite vol- 
ume bias the camera chip needs to be extremely big, 
which is currently a technical obstacle. Alternative se- 
tups where different camera systems focus on the small 
and the large scales simultaneous could be the solution 
to this problem. DNS does not suffer from finite volume 
bias. On the other hand DNS may still have problems 
with small time lags due to interpolation from the Eule- 
rian flow field [H. 

In particular we have studied one-particle statistics in 
terms of Lagrangian structure functions. We have looked 
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at the small time scale behavior which seems to be af- 
fected by the large-scale inhomogeneities present in our 
flow. This led us to work with isotropic ensembles. In 
this way we can test our data against theories developed 
in an isotropic frame such as the multifractal model. 

We do not observe any signs of an inertial range and 
K41 scaling but by Extended Self-Similarity we are able 
to extract a quantitative measure of the structure func- 
tions of high order. From the local slopes of these we cal- 
culate the Lagrangian anomalous scaling exponents and 
find excellent agreement with already published results. 

Measured local slopes of the Lagrangian structure 
functions are quite similar to results obtained with the 
multifractal model for p = 4. With the assumptions and 
physical reasoning leading to the development of the mul- 
tifractal model in mind, this is actually a bit surprising. 
Many of the crucial assumptions behind the multifractal 
model in the Lagrangian frame are not fulfilled. First of 
all the multifractal model is motivated by the invariance 
of the Navier-Stokes equation to an infinite number of 
scaling groups in the limit of infinite Reynolds number 
far from present in our data. Second, an exact result 
such as the 4/5 law does not exist in the Lagrangian 
frame. The phenomenological picture of the multifractal 
model with the flow region consisting of active and inac- 
tive regions in direct contrast to the Richardson picture 



where eddies are space filling does, however, fit observed 
flow features such as a dip region in (r) and anomalous 
scaling. 

For C6(''') and Cs{t) the situation is different. The mul- 
tifractal model fails to describe the data: although the 
qualitative behavior is similar to C4(''') fti^ shift in dip 
position as a function of p is much more pronounced in 
our data than in the multifractal model. We are curious 
to see results from other experiments and DNS simula- 
tions for p > 4. The simple bridging relation, where a 
time lag is solely associated to the time scale, r for the 
local eddy of size r, assumed in the multifractal model 
might be too crude. Also the filtering was shown to have 
an effect on the dip region: the strength for all orders of 
p while the position was only dependent on filtering up 
to p = 6. Since almost all data sets, both from physical 
experiments and from DNS, are noisy and hence must be 
filtered, conclusions on Cpi''') should therefore be made 
with great care. 
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